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Abstract 

Following the monodromy technique performed by Motl and Neitzke, we consider the analytic 
determination of the highly damped (asymptotic) quasi-normal modes of small Schwarzschild-de 
Sitter (SdS) black holes. We comment the result as compared to the recent numerical data of 
Konoplya and Zhidenko. 
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For intermediate times, perturbations of a black hole are dominated by characteristic, 
damped oscillations - the so-called quasi-normal modes (QNMs), the pure tons it emitts in 
the final stage of its perturbation (for reviews, see 

The search for large imaginary, highly damped QNM frequencies has recently attracted a 
great deal of interest. This has been motivated by the work of Hod, who, based on the Bohr 
correspondence principle, proposed that the real part of the asymptotic, highly damped 
QNM frequencies of a black hole can be used to determine the spacing (fundamental unit) 
of its quantum area spectrum |2j. The quantization of black holes horizon area has been 
proposed by Bekenstein thirty years ago based on the adiabatic invariance of the horizon 
area and the analogy with the quantum mechanics of adiabatic systems (for a review, see 
^J). In 3|, Hod has used the highly damped spectrum of QNMs of a Schwarzschild black 
hole, numerically obtained by Nollert e 5] and later confirmed by Andersson [6] and Liu Q. 



Recently, Dreyer |8[ used Hod's conjecture in the context of Loop Quantum Gravity to 
fix a free-parameter (the Barbero-Immirzi parameter) which appears in the area spectrum 
provided by this theory (for reviews, see [2|). In this way, he could naturally obtain the 
Bekenstein-Hawking entropy (see also [10]). Thereafter, several works concerning the highly 



damped QNMs and their role in the black hole area quantization have appeared 



review and more references, see 



1]| (for a 



12|). There have also been studies along the same lines, 



but concerning non-asymptotically flat black holes, namely, black holes in asymptotically 
de Sitter (dS) and anti-de Sitter (AdS) space-times Furthermore, the study of QNMs 
of dS [3| and AdS [3] black holes has been under intensive study. Finally, we mention 
that the works of Motl and Motl and Neitzke ^1 deserve special attention, since they 
first provided an analytic determination of the large imaginary QNMs of a Schwarzschild 
black hole, in agreement with the numerical formula of Nollert la]. Recently, following the 
approach of Motl and Neitzke in Krasnov and Solodukhin 3] rederived the real part of 
the highly damped QNMs of a Schwarzschild black hole by means of a more detailed study, 
on the complex plane, of the relationship between the monodromy associated to the type of 
the differential equation which appears in black hole perturbation problem and the related 
Riemann surfaces. 

The Schwarzschild-de Sitter (SdS) black hole metric in Schwarzschild coordinates 

(we will be using natural units, G — c — k — h— 1, throughout 
the text) ds 2 = -f{r)dt 2 + f(r)- l dr 2 + r 2 {dB 2 + sin 2 dd<p 2 ) , where f{r) = 1 — — — ^ . M 
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is the black hole mass and a the "cosmological radius", related to the cosmological constant 
A>0bya 2 = 3/A. Describing a black hole in a cosmological background, M has to be 
limited to the interval < M < M N = a/3y/3. In this interval, the function f(r) has three 
distinct real zeros : at r = r e , the black hole event horizon, r c ( > r e ) , the cosmological 
event horizon, and tq = — (r e + r c ) . The extreme value Mn gives the Nariai limit [20], for 
which the event and cosmological horizons coincide. 

Associated to each horizon rj there is a surface gravity Ki , defined by Ki = i|^| r=r . . 
Explicitly, we have n e = (r c - r e )(r e - r )/2a 2 r e , k c = (r c - r e )(r c - r )/2a 2 r c , and k = 
(?"e — fo){ r c — ^o)/2a 2 (— ro) ■ The metric parameters M and a can be written in terms of r e 
and r c as a 2 = r 2 + r 2 + r e r c and 2Ma 2 = r e r c {r e + r c ) . 

Assuming a time dependence e tujt , expanding the perturbating field into multipoles and 
defining the tortoise coordinate x by dx/dr = l/f(r) , one can show that the equations for 
scalar, gravitational and electromagnetic perturbations of a SdS black hole are summarized 
by H 

+ ui 2 - V(r(x)) ip(r) = , (1) 



d 2 ip(r) 



dx 2 

where the potential V(r) is 

t + 1) 2M(j 



V(r) = f(r) 



ty 2 ^3 



(2) 



for gravitational (j = 2) and electromagnetic (j = 1) perturbations. £ and j are the 
multipole index and the spin of the perturbating field, respectively. We will consider here 
only electromagnetic and axial gravitational perturbations (in four dimensions, axial and 
polar gravitational perturbations are isospectral). 
Explicitly, the tortoise coordinate is given by 

^J-lnfl-lUJ-lnfl-ZlU-Llnfl-^ . (3) 



2n e \r e J 2k c \ r c J 2kq \ Tq, 

Since V(x) = V(r(x)) goes to zero exponentially for x — * ±oo , in terms of the tortoise 
coordinate the boundary conditions for the QNMs of a SdS black hole are typically the same 
as those of an asymptotically Schwarzschild black hole, namely, ingoing waves at the black 
hole horizon (x — ► — oo) and outgoing waves at the cosmological horizon (x — ► +oo), i.e., in 
view of the dependence e lujt we have 

V>~e^ x ' , x^±oo . (4) 
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Motl and Neitzke |17( have considered the analytic extension to the complex r-plane of 
the perturbation equation (similar to for the Schwarzschild metric. Since in this case 
the equation has (regular) singular points at the origin (r = 0) and at the black hole horizon 
(r = r e ), its solution ip(r) is multivalued around these points. Such multivaluedness is a 
crucial feature in the analysis of 171 ]. for the key idea is to use the QNM boundary conditions 
to compute the monodromy of ip(r) as we run along a conveniently chosen contour in the 
r-plane. Comparison of the monodromy calculated locally (by solving the perturbation 
equation directly) and globally (by direct application of the QNM boundary conditions) 
allowed the determination of the highly damped quasi- normal modes, i.e., those with Imuo > 
> Reu . 

Although the physical region of interest for the problem of QNMs of a SdS black hole is 
r e < r < r c , an analytic extension to the region < r < r c is essential for a "monodromic 
analysis". As we will verify later, the Schwarzschild and SdS cases have similar features 
which concern the behaviour of the potential near the black hole singularity (r = ). 

Due to the logarithmic terms in (JHJ), the analytic continuation to complex-r implies that 
x is a multivalued function of r . But Re x is free of ambiguity and we can determine its sign 
in the r-plane. The multivaluedness of x(r) is also fundamental in this approach. Together 
with the multivaluedness of ip(r) around the (regular) singular point r e , it will lead us to 
the results we look for. 

As in [12I, we also conveniently introduce the variable z = x + const. Choosing the 
constant such that z = for r = , and the branch n = for ln(±l) , from (j3J) we have 
z = x — ni/2K e . 

By analytically continuing r we can distinguish the functions e T%U)X by their monodromy 
at r = r e . This enables us to define the boundary condition at r = r e by simply requiring 
that ip(r) have monodromy e 2lt i u l 2K ^) — e ™/K e on a clockwise contour around this point. 
We note that we cannot do the same at r = r c , since as r is now complex, x does not 
tend to +00 for r = r c . Nevertheless, if we consider the special case of small black holes 
(r c >> r e ), we can approximately consider that x >> 1 for r = r c , as an expansion of x(r) 
for (r c /r e ) » 1 shows. 

n 

Now we can closely follow [17J] and in order to define the boundary condition at r c , we 
simply analytically continue ip{r) via rotation to the line Im (ujx) = . Hence we will 
have a purely oscillatory asymptotic behaviour on the real line ujx , what implies that it 
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is possible to select a particular solution by specifying its asymptotics. We will match the 
asymptotics along the line Re x = (corresponding to the left part of the contour 7 between 
A and B in the figure and use the boundary conditions to determine the QNMs. Taking 
the limit where uj is almost purely imaginary (highly damped modes), we see that the line 
Im {ujx) = is slightly sloped off the line Re x = . Since there are two possible directions 
for the rotation, we can choose the one corresponding to an angle smaller than 7r/2 . 

From (j3J), we note that for QNMs we have Imu > . Assuming initially that Re to > , 
in view of the condition Im (ujx) = , we see that x = +00 is rotated to ujx = +00 , and on 
this line the boundary condition at x = +00 is 

V>(r) ~ e~ iuJX , ujx -> +00 . (5) 

Let us first compute the monodromy locally. This will be done by matching the asymp- 
totics along the line Im (ujx) = (or along Im (ujz) = , for uj (almost) purely imaginary). 
We will start at A and move along the contour 7 towards the origin (see figure Q). The 
asymptotics at A can be matched to that near the origin. First, to obtain the solution near 
the origin, we expand x near r = and thus obtain for small z 

z^ f3r 2 , (6) 

where p=\(-L?--L^- . 

Therefore, near the origin the potential behaves as 

V(z)k — ^ ^ . (7) 

We remark that © has been deduced without using the approximation of small black 
holes. Hence, the near origin behaviour (J2j) remains valid in general. 

In view of (0), (JTJ can be reduced to a Bessel equation. Its solution near the origin is 
then given by 

ip(z) = A + c + y/uJzJ u (ujz) + A^c^\fujz.I- v (ujz) , (8) 
with the index v being given by 

u = ^l + (4M(3y(f-l) ■ (9) 

We will now match the solution (jHJ) to the asymptotics away from the origin. Taking the 
asymptotic behaviour of J± u (uuz) as ujz — > 00 , we can choose the constants c± in (jHJ) such 



FIG. 1: The complex r-plane and the contour 7. The hachured part marks the region where 
Rex < 0. 



that 

c±\fujzj± v {u)z) ~ 2 cos(W — a±) as uz — > 00 , (10) 

where a± = |(1 ± 2v) . 

From (|T0|) and (jHJ) , and making use of the boundary condition (JHJ) , we can write for the 
asymptotics at A 

?P{z) ~ {A + e ia+ + A_e ia -)e~ luJZ as ujz -> 00 , (11) 

since 

A + e^ Q + + A_e" iQ - = . (12) 

In order to follow the contour until i? we have to turn around the origin. Thus we 
must perform a rotation of 3tc/2 in the r-plane, corresponding to 3tc in the z-plane. Using 
the asymptotic behaviour of the Bessel function near the origin, we can match it to the 
asymptotic solution. In fact, using for the Bessel function the behaviour J± v (z) ~ (|) 
near the origin, after the 3n rotation we have c±y/uJzJ± u (tuz) ~ e 6ta± 2 cos(— cuz — a±) , as 
ujz — > — 00 , which combined with (jHJ) gives for the aymptotics at B 

il){z) ~ (A + e 5ia+ + A_e 5ia -)e- iu)Z + (A + e 7ia+ + A_e 7ia -)e iuz as uz -> -00 . (13) 
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We can go from B back to A along the large semi-circle (|r| ~ \r c \)- In this region, taking 
into account our approximation of small black holes, we have x » 1 , and thus V(x) goes to 
zero, such that u 2 dominates the potential, and we can approximate the solution of the wave 
equation there by plane waves. Hence, the coefficient of e~ lu}Z remains essentially unchanged 
as we return to A , whereas the same cannot be said for the coefficient of e luJZ , which gives 
an exponentially small contribution to ip(r) in the region where Rex > . Therefore, from 
(fTTj) and (fH?|) . we see that the monodromy around the contour 7 is simply given by 

A 5ia+ + A_e 5ia ~ . 

— — - — ; = -(1 + 2 cos 2m/) , 14 

A + e ia + + A_e ia - 

where we have taken into account ()12|1 . 

The global computation of the monodromy around 7 can be obtained by making use of 

the QNM boundary condition at r e , since the only singularity of ip(r) or e~ lu3Z inside the 

contour 7 occurs at this point. The boundary condition at r e implies that after a full trip 

around the contour clockwise, ip{r) acquires a phase e nul ^ Ke , while e~ wz acquires a phase 

e -TTUj/K e _ Therefore, we must multiply the coefficient of e~ tuJZ in the asymptotics of tp(r) by 

Comparing the local monodromy (fT4*|) with the global one, we find 

e 2 ™"/** = -(1 + 2 cos 2m/) . (15) 

For the choice Re u < we would have to consider the rotation in the opposite direction 
and then reverse considerations on the points A and B , as well as to run the contour 7 
in the counter-clockwise direction. We thus conclude that the asymptotic, highly damped 
quasi-normal mode spectrum of small SdS black holes is given by 

= ( n + -)i ± — In |1 + 2cos2vrz/| , (16) 
n e 2 2n 

with n — > 00 . 

From the relation AM 3 = — 1 + n , re J~ c , 2 , if we take the limit as r r — ► 00 , which 
corresponds to A = , we find AM 3 = — 1 , such that we recover the result of Motl and 
Neitzke for asymptotically flat Schwarzschild black holes [171 ]. 

8irMuj n = (2n + l)m ± In |1 + 2cos7rj| . (17) 

The result ()16|) shows that the real part of the asymptotic modes has a constant value 
and then it fails to reproduce the oscillatory behaviour recently found by Konoplya and 



Zhidenko 



22j in the case of gravitational perturbations (j = 2). Nevertheless, for the 



electromagnetic case (j = 1 ), Reu n — , what agrees with the result of those authors, but 
the behaviour of the imaginary part of the modes, as given by only partially agrees 
with that found in 22] , since it does not show an oscillatory behaviour in its spacing, found 
for both gravitational and electromagnetic cases in (22}]. Accorging to [2^, an oscillatory 
term should appear in Imu n , in addition to (n + \)n e - 

Summarizing, by applying the monodromy technique as performed by Motl and Neitzke 
^|, we have considered the determination of the asymptotic, highly damped quasi-normal 
mode spectrum of small (r e << r c ) Schwarzschild-de Sitter black holes. However, we note 
that before we have restricted the analysis for small black holes, we verified that the dominant 
contribution for potential near the black hole singularity is exactly of the same type as found 
for the asymptotically flat case, i.e., the solution is a Bessel function around the origin. This 
result is important in applying the monodromy approach to find the asymptotic QNMs for 
general SdS black holes. Compared to the numerical results of Konoplya and Zhidenko [22] , 
it seems our results do not give the correct asymptotic QNM frequencies for small SdS black 
holes. Consequently, in the context of the monodromy technique, it is important to search 



22 



for the general analytic solution and compare it with the available numerical data in 
n fact, this has already been done very recently by Cardoso et al. and it will appear soon 
2J. 
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